Physical Phenomena Governing Mineral Morphogenesis in Molluscan Nacre.

Mollusks, as well as many other living organisms, have the ability to shape mineral crystals into unconventional morphologies and to assemble them into complex functional mineral-organic structures, an observation that inspired tremendous research efforts in scientific and technological domains. Despite these, a biochemical toolkit that accounts for the formation of the vast variety of the observed mineral morphologies cannot be identified yet. Herein, phase-field modeling of molluscan nacre formation, an intensively studied biomineralization process, is used to identify key physical parameters that govern mineral morphogenesis. Manipulating such parameters, various nacre properties ranging from the morphology of a single mineral building block to that of the entire nacreous assembly are reproduced. The results support the hypothesis that the control over mineral morphogenesis in mineralized tissues happens via regulating the physico-chemical environment, in which biomineralization occurs: the organic content manipulates the geometric and thermodynamic boundary conditions, which in turn, determine the process of growth and the form of the biomineral phase. The approach developed here has the potential of providing explicit guidelines for the morphogenetic control of synthetically formed composite materials.

Phase-Field Modeling of Biomineralization in Mollusks and Corals: Microstructure vs Formation Mechanism.

While biological crystallization processes have been studied on the microscale extensively, there is a general lack of models addressing the mesoscale aspects of such phenomena. In this work, we investigate whether the phase-field theory developed in materials' science for describing complex polycrystalline structures on the mesoscale can be meaningfully adapted to model crystallization in biological systems. We demonstrate the abilities of the phase-field technique by modeling a range of microstructures observed in mollusk shells and coral skeletons, including granular, prismatic, sheet/columnar nacre, and sprinkled spherulitic structures. We also compare two possible micromechanisms of calcification: the classical route, via ion-by-ion addition from a fluid state, and a nonclassical route, crystallization of an amorphous precursor deposited at the solidification front. We show that with an appropriate choice of the model parameters, microstructures similar to those found in biomineralized systems can be obtained along both routes, though the time-scale of the nonclassical route appears to be more realistic. The resemblance of the simulated and natural biominerals suggests that, underneath the immense biological complexity observed in living organisms, the underlying design principles for biological structures may be understood with simple math and simulated by phase-field theory.

Crystal nucleation and growth of spherulites demonstrated by coral skeletons and phase-field simulations.

Spherulites are radial distributions of acicular crystals, common in biogenic, geologic, and synthetic systems, yet exactly how spherulitic crystals nucleate and grow is still poorly understood. To investigate these processes in more detail, we chose scleractinian corals as a model system, because they are well known to form their skeletons from aragonite (CaCO3) spherulites, and because a comparative study of crystal structures across coral species has not been performed previously. We observed that all 12 diverse coral species analyzed here exhibit plumose spherulites in their skeletons, with well-defined centers of calcification (CoCs), and crystalline fibers radiating from them. In 7 of the 12 species, we observed a skeletal structural motif not observed previously: randomly oriented, equant crystals, which we termed "sprinkles". In Acropora pharaonis, these sprinkles are localized at the CoCs, while in 6 other species, sprinkles are either layered at the growth front (GF) of the spherulites, or randomly distributed. At the nano- and micro-scale, coral skeletons fill space as much as single crystals of aragonite. Based on these observations, we tentatively propose a spherulite formation mechanism in which growth front nucleation (GFN) of randomly oriented sprinkles, competition for space, and coarsening produce spherulites, rather than the previously assumed slightly misoriented nucleations termed "non-crystallographic branching". Phase-field simulations support this mechanism, and, using a minimal set of thermodynamic parameters, are able to reproduce all of the microstructural variation observed experimentally in all of the investigated coral skeletons. Beyond coral skeletons, other spherulitic systems, from aspirin to semicrystalline polymers and chocolate, may also form according to the mechanism for spherulite formation proposed here. STATEMENT OF SIGNIFICANCE: Understanding the fundamental mechanisms of spherulite nucleation and growth has broad ranging applications in the fields of metallurgy, polymers, food science, and pharmaceutical production. Using the skeletons of reef-building corals as a model system for investigating these processes, we propose a new spherulite growth mechanism that can not only explain the micro-structural diversity observed in distantly related coral species, but may point to a universal growth mechanism in a wide range of biologically and technologically relevant spherulitic materials systems.

Crystal growth kinetics as an architectural constraint on the evolution of molluscan shells.

Molluscan shells are a classic model system to study formation-structure-function relationships in biological materials and the process of biomineralized tissue morphogenesis. Typically, each shell consists of a number of highly mineralized ultrastructures, each characterized by a specific 3D mineral-organic architecture. Surprisingly, in some cases, despite the lack of a mutual biochemical toolkit for biomineralization or evidence of homology, shells from different independently evolved species contain similar ultrastructural motifs. In the present study, using a recently developed physical framework, which is based on an analogy to the process of directional solidification and simulated by phase-field modeling, we compare the process of ultrastructural morphogenesis of shells from 3 major molluscan classes: A bivalve Unio pictorum, a cephalopod Nautilus pompilius, and a gastropod Haliotis asinina We demonstrate that the fabrication of these tissues is guided by the organisms by regulating the chemical and physical boundary conditions that control the growth kinetics of the mineral phase. This biomineralization concept is postulated to act as an architectural constraint on the evolution of molluscan shells by defining a morphospace of possible shell ultrastructures that is bounded by the thermodynamics and kinetics of crystal growth.

Biomineralization as a Paradigm of Directional Solidification: A Physical Model for Molluscan Shell Ultrastructural Morphogenesis.

Molluscan shells are a model system to understand the fundamental principles of mineral formation by living organisms. The diversity of unconventional mineral morphologies and 3D mineral-organic architectures that comprise these tissues, in combination with their exceptional mechanical efficiency, offers a unique platform to study the formation-structure-function relationship in a biomineralized system. However, so far, morphogenesis of these ultrastructures is poorly understood. Here, a comprehensive physical model, based on the concept of directional solidification, is developed to describe molluscan shell biomineralization. The capacity of the model to define the forces and thermodynamic constraints that guide the morphogenesis of the entire shell construct-the prismatic and nacreous ultrastructures and their transitions-and govern the evolution of the constituent mineralized assemblies on the ultrastructural and nanostructural levels is demonstrated using the shell of the bivalve Unio pictorum. Thereby, explicit tools for novel bioinspired and biomimetic bottom-up materials design are provided.

Grain coarsening in two-dimensional phase-field models with an orientation field.

In the literature, contradictory results have been published regarding the form of the limiting (long-time) grain size distribution (LGSD) that characterizes the late stage grain coarsening in two-dimensional and quasi-two-dimensional polycrystalline systems. While experiments and the phase-field crystal (PFC) model (a simple dynamical density functional theory) indicate a log-normal distribution, other works including theoretical studies based on conventional phase-field simulations that rely on coarse grained fields, like the multi-phase-field (MPF) and orientation field (OF) models, yield significantly different distributions. In a recent work, we have shown that the coarse grained phase-field models (whether MPF or OF) yield very similar limiting size distributions that seem to differ from the theoretical predictions. Herein, we revisit this problem, and demonstrate in the case of OF models [R. Kobayashi, J. A. Warren, and W. C. Carter, Physica D 140, 141 (2000)PDNPDT0167-278910.1016/S0167-2789(00)00023-3; H. Henry, J. Mellenthin, and M. Plapp, Phys. Rev. B 86, 054117 (2012)PRBMDO1098-012110.1103/PhysRevB.86.054117] that an insufficient resolution of the small angle grain boundaries leads to a log-normal distribution close to those seen in the experiments and the molecular scale PFC simulations. Our paper indicates, furthermore, that the LGSD is critically sensitive to the details of the evaluation process, and raises the possibility that the differences among the LGSD results from different sources may originate from differences in the detection of small angle grain boundaries.

Topological defects in two-dimensional orientation-field models for grain growth.

Standard two-dimensional orientation-field-based phase-field models rely on a continuous scalar field to represent crystallographic orientation. The corresponding order parameter space is the unit circle, which is not simply connected. This topological property has important consequences for the resulting multigrain structures: (i) trijunctions may be singular; (ii) for each pair of grains there exist two different grain boundary solutions that cannot continuously transform to one another; (iii) if both solutions appear along a grain boundary, a topologically stable, singular point defect must exist between them. While (i) can be interpreted in the classical picture of grain boundaries, (ii) and therefore (iii) cannot. In addition, singularities cause difficulties, such as lattice pinning in numerical simulations. To overcome these problems, we propose two formulations of the model. The first is based on a three-component unit vector field, while in the second we utilize a two-component vector field with an additional potential. In both cases, the additional degree of freedom introduced makes the order parameter space simply connected, which removes the topological stability of these defects.

Ternary eutectic dendrites: Pattern formation and scaling properties.

Extending previous work [Pusztai et al., Phys. Rev. E 87, 032401 (2013)], we have studied the formation of eutectic dendrites in a model ternary system within the framework of the phase-field theory. We have mapped out the domain in which two-phase dendritic structures grow. With increasing pulling velocity, the following sequence of growth morphologies is observed: flat front lamellae → eutectic colonies → eutectic dendrites → dendrites with target pattern → partitionless dendrites → partitionless flat front. We confirm that the two-phase and one-phase dendrites have similar forms and display a similar scaling of the dendrite tip radius with the interface free energy. It is also found that the possible eutectic patterns include the target pattern, and single- and multiarm spirals, of which the thermal fluctuations choose. The most probable number of spiral arms increases with increasing tip radius and with decreasing kinetic anisotropy. Our numerical simulations confirm that in agreement with the assumptions of a recent analysis of two-phase dendrites [Akamatsu et al., Phys. Rev. Lett. 112, 105502 (2014)], the Jackson-Hunt scaling of the eutectic wavelength with pulling velocity is obeyed in the parameter domain explored, and that the natural eutectic wavelength is proportional to the tip radius of the two-phase dendrites. Finally, we find that it is very difficult/virtually impossible to form spiraling two-phase dendrites without anisotropy, an observation that seems to contradict the expectations of Akamatsu et al. Yet, it cannot be excluded that in isotropic systems, two-phase dendrites are rare events difficult to observe in simulations.

Heterogeneous nucleation of/on nanoparticles: a density functional study using the phase-field crystal model.

Crystallization of supersaturated liquids usually starts by heterogeneous nucleation. Mounting evidence shows that even homogeneous nucleation in simple liquids takes place in two steps; first a dense amorphous precursor forms, and the crystalline phase appears via heterogeneous nucleation in/on the precursor cluster. Herein, we review recent results by a simple dynamical density functional theory, the phase-field crystal model, for (precursor-mediated) homogeneous and heterogeneous nucleation of nanocrystals. It will be shown that the mismatch between the lattice constants of the nucleating crystal and the substrate plays a decisive role in determining the contact angle and nucleation barrier, which were found to be non-monotonic functions of the lattice mismatch. Time dependent studies are essential as investigations based on equilibrium properties often cannot identify the preferred nucleation pathways. Modeling of these phenomena is essential for designing materials on the basis of controlled nucleation and/or nano-patterning.

Heterogeneous crystal nucleation: the effect of lattice mismatch.

A simple dynamical density functional theory is used to investigate freezing of an undercooled liquid in the presence of a crystalline substrate. We find that the adsorption of the crystalline phase on the substrate, the contact angle, and the height of the nucleation barrier are nonmonotonic functions of the lattice constant of the substrate. We show that the free-growth-limited model of particle-induced freezing by Greer et al. [Acta Mater. 48, 2823 (2000)] is valid for larger nanoparticles and a small anisotropy of the interface free energy. Faceting due to the small size of the foreign particle or a high anisotropy decouples free growth from the critical size of homogeneous nuclei.

Amorphous nucleation precursor in highly nonequilibrium fluids.

Dynamical density-functional simulations reveal structural aspects of crystal nucleation in undercooled liquids: The first appearing solid is amorphous, which promotes the nucleation of bcc crystals but suppresses the appearance of the fcc and hcp phases. These findings are associated with features of the effective interaction potential deduced from the amorphous structure.

Polymorphism, crystal nucleation and growth in the phase-field crystal model in 2D and 3D.

We apply a simple dynamical density functional theory, the phase-field crystal (PFC) model of overdamped conservative dynamics, to address polymorphism, crystal nucleation, and crystal growth in the diffusion-controlled limit. We refine the phase diagram for 3D, and determine the line free energy in 2D and the height of the nucleation barrier in 2D and 3D for homogeneous and heterogeneous nucleation by solving the respective Euler-Lagrange (EL) equations. We demonstrate that, in the PFC model, the body-centered cubic (bcc), the face-centered cubic (fcc), and the hexagonal close-packed structures (hcp) compete, while the simple cubic structure is unstable, and that phase preference can be tuned by changing the model parameters: close to the critical point the bcc structure is stable, while far from the critical point the fcc prevails, with an hcp stability domain in between. We note that with increasing distance from the critical point the equilibrium shapes vary from the sphere to specific faceted shapes: rhombic dodecahedron (bcc), truncated octahedron (fcc), and hexagonal prism (hcp). Solving the equation of motion of the PFC model supplied with conserved noise, solidification starts with the nucleation of an amorphous precursor phase, into which the stable crystalline phase nucleates. The growth rate is found to be time dependent and anisotropic; this anisotropy depends on the driving force. We show that due to the diffusion-controlled growth mechanism, which is especially relevant for crystal aggregation in colloidal systems, dendritic growth structures evolve in large-scale isothermal single-component PFC simulations. An oscillatory effective pair potential resembling those for model glass formers has been evaluated from structural data of the amorphous phase obtained by instantaneous quenching. Finally, we present results for eutectic solidification in a binary PFC model.

Phase field theory of heterogeneous crystal nucleation.

The phase field approach is used to model heterogeneous crystal nucleation in an undercooled pure liquid in contact with a foreign wall. We discuss various choices for the boundary condition at the wall and determine the properties of critical nuclei, including their free energy of formation and the contact angle as a function of undercooling. For particular choices of boundary conditions, we may realize either an analog of the classical spherical cap model or decidedly nonclassical behavior, where the contact angle decreases from its value taken at the melting point towards complete wetting at a critical undercooling, an analogue of the surface spinodal of liquid-wall interfaces.

Multiscale approach to CO2 hydrate formation in aqueous solution: phase field theory and molecular dynamics. Nucleation and growth.

A phase field theory with model parameters evaluated from atomistic simulations/experiments is applied to predict the nucleation and growth rates of solid CO(2) hydrate in aqueous solutions under conditions typical to underwater natural gas hydrate reservoirs. It is shown that under practical conditions a homogeneous nucleation of the hydrate phase can be ruled out. The growth rate of CO(2) hydrate dendrites has been determined from phase field simulations as a function of composition while using a physical interface thickness (0.85+/-0.07 nm) evaluated from molecular dynamics simulations. The growth rate extrapolated to realistic supersaturations is about three orders of magnitude larger than the respective experimental observation. A possible origin of the discrepancy is discussed. It is suggested that a kinetic barrier reflecting the difficulties in building the complex crystal structure is the most probable source of the deviations.

Growth and form of spherulites.

Many structural materials (metal alloys, polymers, minerals, etc.) are formed by quenching liquids into crystalline solids. This highly nonequilibrium process often leads to polycrystalline growth patterns that are broadly termed "spherulites" because of their large-scale average spherical shape. Despite the prevalence and practical importance of spherulite formation, only rather qualitative concepts of this phenomenon exist. It is established that phase field methods naturally account for diffusional instabilities that are responsible for dendritic single-crystal growth. However, a generalization of this model is required to describe spherulitic growth patterns, and in the present paper we propose a minimal model of this fundamental crystal growth process. Our calculations indicate that the diversity of spherulitic growth morphologies arises from a competition between the ordering effect of discrete local crystallographic symmetries and the randomization of the local crystallographic orientation that accompanies crystal grain nucleation at the growth front [growth front nucleation (GFN)]. This randomization in the orientation accounts for the isotropy of spherulitic growth at large length scales and long times. In practice, many mechanisms can give rise to GFN, and the present work describes and explores three physically prevalent sources of disorder that lead to this kind of growth. While previous phase field modeling elucidated two of these mechanisms--disorder created by particulate impurities or other static disorder or by the dynamic heterogeneities that spontaneously form in supercooled liquids (even pure ones)--the present paper considers an additional mechanism, crystalline branching induced by a misorientation-dependent grain boundary energy, which can significantly affect spherulite morphology. We find the entire range of observed spherulite morphologies can be reproduced by this generalized phase field model of polycrystalline growth.

A general mechanism of polycrystalline growth.

Most research into microstructure formation during solidification has focused on single-crystal growth ranging from faceted crystals to symmetric dendrites. However, these growth forms can be perturbed by heterogeneities, yielding a rich variety of polycrystalline growth patterns. Phase-field simulations show that the presence of particulates (for example, dirt) or a small rotational-translational mobility ratio (characteristic of high supercooling) in crystallizing fluids give rise to similar growth patterns, implying a duality in the growth process in these structurally heterogeneous fluids. Similar crystallization patterns are also found in thin polymer films with particulate additives and pure films with high supercooling. This duality between the static and dynamic heterogeneity explains the ubiquity of polycrystalline growth patterns in polymeric and other complex fluids.

Growth of 'dizzy dendrites' in a random field of foreign particles.

Microstructure plays an essential role in determining the properties of crystalline materials. A widely used method to influence microstructure is the addition of nucleating agents. Observations on films formed from clay-polymer blends indicate that particulate additives, in addition to serving as nucleating agents, may also perturb crystal growth, leading to the formation of irregular dendritic morphologies. Here we describe the formation of these 'dizzy dendrites' using a phase-field theory, in which randomly distributed foreign particle inclusions perturb the crystallization by deflecting the tips of the growing dendrite arms. This mechanism of crystallization, which is verified experimentally, leads to a polycrystalline structure dependent on particle configuration and orientation. Using computer simulations we demonstrate that additives of controlled crystal orientation should allow for a substantial manipulation of the crystallization morphology.

Nucleation and bulk crystallization in binary phase field theory.

We present a phase field theory for binary crystal nucleation. In the one-component limit, quantitative agreement is achieved with computer simulations (Lennard-Jones system) and experiments (ice-water system) using model parameters evaluated from the free energy and thickness of the interface. The critical undercoolings predicted for Cu-Ni alloys accord with the measurements, and indicate homogeneous nucleation. The Kolmogorov exponents deduced for dendritic solidification and for "soft impingement" of particles via diffusion fields are consistent with experiment.